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abstract 

We analyze when and why unitarity violations might occur in quantum cosmology 
restricted to minisuperspace. To this end we discuss in detail backscattering transitions 
between expanding and contracting solutions of the Wheeler-DeWitt equation. We first 
show that upon neglecting only backscattering, one obtains an intermediate regime in 
which matter evolves unitarily but which does not correspond to any Schrodinger equa- 
tion in a given geometry since gravitational backreaction effects are taken into account 
at the quantum level. We then show that backscattering amplitudes are exponentially 
smaller than matter transition amplitudes. Both results follow from an adiabatic treat- 
ment valid for macroscopic universes. To understand how backscattering and the inter- 
mediate regime should be interpreted, we review the problem of electronic transitions 
induced by nuclear motion since it is mathematically very similar. In this problem, 
transition amplitudes are obtained from the conserved current. The same applies to 
quantum cosmology and indicates that probability amplitudes should based on the cur- 
rent when backscattering is neglected. We then review why, in a relativistic context, 
backscattering is interpreted as pair production whereas it is not in the non relativistic 
case. In each example the correct interpretation is obtained by coupling the system 
to an external quantum device. From the absence of such external systems in cosmol- 
ogy, we conclude that backscattering does not have a unique consistent interpretation 
in quantum cosmology. 
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1 Introduction 



Quantum gravity suggests the possibility of non unitary evolution. In particular singu- 
larities and regions of strong curvature, such as black holes and the Big Bang, are natural 
places where such effects could occur. However, in these cases, the lack of comprehension 
of Planck scale physics precludes us from reaching any definitive conclusions. 

More surprisingly, when the universe is macroscopic and the curvature is small, the 
possibility of non unitarity remains. Moreover, in this long wavelength regime, it seems 
unavoidable that quantum gravity implies the Wheeler-DeWitt (WDW) equation]]] 
(truncated to a finite set of modes so as to eliminate ultraviolet problems) since this 
constraint expresses and guarantees reparametrisation invariance. 

Thus the theoretical basis for investigating unitarity violation in quantum cosmology 
seems well established. Nevertheless no definite analysis has been carried out, even in 
minisuperspace. Indeed, in the literature one finds either prudent discourses which avoid 
the problem or claims ranging from "there are no unitary violations" || to "there are 
finite violations of unitarity" || . 

In order to understand the origin of the difficulties which have hindered previous work 
on this question, it is necessary to recall how the usual unitary evolution of quantum 
matter in a given four-geometry is recovered from the WDW equation. This has already 
been thoroughly investigated P]- [p^] and it is now well understood that if the wave 



function is in a tight wave packet, i.e. the spread of matter energy is small, and if 
gravity can be treated semiclassicaly, i.e. it can be described by a WKB wave function, 
then one should factor out this WKB wave and interpret the remainder as the wave 
function of matter which evolves according to the Schrodinger equation. 

Therefore, when investigating the problem of unitary violations, two problems must 
be confronted. The first one is technical and concerns the precise mathematical char- 
acterization of the corrections to the Schrodinger equation which are due to the aban- 
donment of the two restrictions. In this respect, it should be mentioned that most 
approaches have kept as starting point (i.e. as an ansatz) that the total wave func- 
tion can be factorized into a gravitational wave characterizing the background and the 
rest which is still interpreted as the matter wave function. It is then very difficult to 
distinguish the intrinsic corrections from the those due to this factorization ansatz. 

The second aspect is more fundamental and has to do with the interpretation of the 
wave function of the universe ^(a, 0), solution of the WDW equation (a designates the 
scale factor and matter variables in a minisuperspace context). In particular, in order 
to address the question of possible violations of unitarity, one must know how to extract 
the wave function of matter when matter is not in a tight wave packet. Several proposals 
have been made in the literature, see M for a review. Two of the most specific ones 
are the conditional probability interpretation wherein the probability to find at a is 
given by |\l/(a,0)| 2 / / e?0|\l/(a, 0) | 2 and the Vilenkin[[|] interpretation which is based on 
the current carried by \l/(a, 0). 

we circum- 



In this paper, by exploiting the generalized adiabatic treatment of [111 



vent the difficulties which have hindered previous work. The major advantage of this 
treatment is to organize the WDW equation in such a way that the corrections to the 
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two approximations can be separately and systematically analyzed before adopting any 
interpretation scheme. We then focus on the possibility of unitary violations in the light 
of the mathematical properties of these corrections. Our analysis will be restricted to 
minisuperspace models. In a subsequent paper, we intend to extend it to cosmologies 
characterized by many variables. 

The first important conclusion of this work (which generalizes that of fl3|| ) is that 
when matter is not in a tight wave packet but in the absence of backscattering between 
waves with positive and negative gravitational momentum (which correspond to con- 
tracting and expanding universes), the only consistent interpretation is that based on 
the current. From our analysis, it will also be clear why when matter is in a tight wave 
packet, i.e. when both approximations mentioned before are exact, different interpreta- 
tions, such as the probability or the current interpretations, become equivalent. 

The second point of this article is to analyze the importance of the effects induced by 
backscattering of gravitational waves. We show that because of destructive interference, 
the amplitude for a forward wave to backscatter is proportional to exp(—d a p/p 2 ) where 
p is the momentum of a. This is much smaller than previous estimates that suggested 
that the amplitudes were proportional to the quantity d a p/p 2 itself. 

The third main conclusion concerns the consequences for matter evolution of the 
couplings induced by backscattering. We show that even though these transition am- 
plitudes can be precisely evaluated, it is impossible to deduce what are exactly their 
physical consequences. That is, quantum cosmology is an incomplete theory that does 
not carry its own interpretation. Happily this incompleteness only manifests itself at 
the level of backscattering and therefore concerns exponentially small effects. 

To reach these conclusions we address successively the mathematical and interpreta- 
tional aspects. In the first part of the article (sections |2] to |5J), we carry out an analysis 
of the WDW equation which does not make the ansatz that matter should be in a tight 



wave packet. Dropping this ansatz has already been advocated in [lC][p 1 J[12|[13|. Here 



we shall follow the treatment of [M (see also the earlier work of [IS]) which does not 
require the validity of the WKB approximation. Thus it can be used to obtain quan- 
titative estimates for both amplitudes to backscatter and effects of widely spread wave 
functions. From this treatment one unambiguously establishes that: 

• When backscattering amplitudes can be neglected and matter is in a tight wave 
packet, matter obeys the Schrodinger equation in the background defined by the 
center of the wave packet, thereby recovering the situation of refs. [B], HI [3J. 

• When matter is no longer in a tight wave packet, but backscattering amplitudes 
are still neglected, matter obeys a Schrodinger-like equation wherein evolution is 
parameterized by a and wherein gravitational back reaction effects are included at 
the quantum level, i.e. not only in the mean. Unitarity follows from the ordered 
nature of the WKB propagation of a. 

• When backscattering amplitudes are not neglected, matter evolution is drastically 
modified in the sense that matter states associated with forward and backward 
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propagating universes interact. This directly follows from the second order char- 
acter of the WDW equation in d a . 

These results follow from the separation of the length scales which govern the dynamics 
in quantum cosmology when the universe is macroscopic. The Compton wave length of 
a, 1/p, is much smaller than both the Compton wave length of a typical matter degree 
of freedom and than the time over which the geometry changes {d a p/p)~ l . 



In the second part of this article (sections |6] to [11]) we address the problem of in- 



terpretation. As basis for the discussion, we consider two physical examples whose 
mathematical description is similar to the WDW equation. These examples are elec- 
tronic non adiabatic transitions during atomic collisions and pair creation in an external 
potential. The manner in which the wave functions should be interpreted in these ex- 
amples is standard text book material. Nevertheless, it is very instructive to delve in 
detail into the mathematical and physical basis of these interpretations in order to see 
if they apply to quantum cosmology. 

In both cases the conserved current serves as a basis for extracting predictions from 
the theory. The basic mathematical reason is that this is the only quantity which 
is conserved by virtue of the equations of motion. In cosmology, in the absence of 
backscattering, this procedure leads unambiguously to the current interpretation. In 
this derivation, the existence of well separated length scales plays also a crucial role: It 
provides the possibility of asking questions concerning light degrees of freedom under the 
condition that the heavy one is propagating forward and located in a small neighborhood. 

The other important point is that when backscattering between the forward and 
backward propagating waves occurs, an additional principle is required in order to inter- 
pret backscattering amplitudes. Moreover, this principle differs for the two examples. In 
each case, it is dictated by the possibility of coupling the system to a quantum external 
device which can act as a measuring device. In quantum cosmology, no such external 
quantum device can be introduced, and it is therefore impossible to decide which of these 
two interpretations, or possibly a completely different third one, applies to the WDW 
equation. Hopefully the resolution of this problem can ultimately be obtained from a 
deeper understanding of the inner (ie. ultraviolet) structure of quantum gravity. 



2 Backscattering 

In view of the importance of backscattering effects in the problem of unitarity violation 
in quantum cosmology, it is appropriate to understand clearly why backscattering occurs 
and how to compute backscattering amplitudes. To this end, we first consider the simple 
problem of a universe in which matter stays in a given eigenstate \n) of constant energy 
E n . This enables us to study the corrections to the WKB approximation unencumbered 
by the details of matter evolution. This analysis will serve as a guide to the next section 
wherein the matter Hamiltonian no longer commutes with the WDW equation. 

With this mathematical simplification, the WDW equation in minisuperspace reduces 
to one degree of freedom moving in a potential 

H T m = [G 2 d 2 a + V{a) + 2aGE n ]m = (1) 
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where G is Newton's constant. The gravitational potential V(a) comes from the cur- 
vature of space (for instance, V = —a 2 + Aa 4 in a closed universe with a cosmological 
constant A). 

The WKB solutions to eq. (jl]) with unit Wronskian are 

e — i J a da'p n (a') 



and its complex conjugate, where p n (a) = ^yV(a) + 2aGE n is the classical solution 
of Ht = when matter has energy E n . These solutions correspond to expanding and 
contracting universes. 

This identification relies on a sign convention which we now explain. First one 
expands (as usual when building wave packets) the phase in eq. @ around a reference 
energy E to obtain 

„—i f a da'p n (a') 

Xn(a)^ e —L=— c -*r-'ft(^ . (3) 



Classical mechanics tells us that ^\e=e = ±4^ where t(a) is the proper time obtained 
from the classical trajectory of a driven by the reference matter energy E. The ± sign 
correspond to expanding and contracting universes respectively. Thus we obtain 

-i J a da'p n (a') 
V ^Pn 

If one adopts the usual convention that the phase dependence on energy and time should 
be e~ tEnt , then one must take the — sign in eq. Hence Xn corresponds to an 

expanding universe and x* n t° a contracting universe. Obviously this choice is purely 
conventional and has nothing intrinsic to it. The problem of defining an intrinsic arrow 
of time in quantum cosmology is a delicate one, and we shall briefly return to it at the 
end of this paper. 

There are two kinds of corrections to the WKB solutions. The first are corrections 
to the phase and norm of Xn and x* n separately. They are obtained by a standard 
calculation in which on writes ^> = Ae iS ^ n , inserts this ansatz into eq. (HD and calculates 



A and S as series in h, see for instance||16||. The second kind of corrections are due 



to the couplings between the forward and the backward solution. They describe the 
backscattering amplitudes engendered by the a dependent effective potential V n (a) = 
V(a) + 2aGE n . 

We are interested in the second corrections since, as we shall later see, these are the 
relevant ones for the problem of unitarity violations. To analyze them we rewrite ^ as 

* = C n (a) X n(a) +V n (a) X * n (a) , (5) 
d a ^ = ip n {a) [C n {a)Xn{o) ~ T> n (a)x* n {a)} ■ (6) 



The coefficients C n and V n are uniquely determined by these equations |L7| and the 
conserved current (Wronskian) takes the simple form 

q*id a y = \C n (a) | 2 - \V n (a) | 2 = constant . (7) 
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Differentiating eq. (|5p and comparing with eq. (|6|) yields 

d a C n Xn{a) + d a V n x* n {a) - —— [C n {a)x n {a) + V n (a)x* n {a)} = . (8) 

2p„ 

Inserting eq. (J6|) into eq. (P and using eq. (||) yields first order coupled equations 

P) r — ^ a P n f a da p n 

2 p n 

daVn = ^^ e +»/-&ft. (9) 

which are equivalent to the original WDW equation ([[]). 

From eq. (||) one reads that the instantaneous coupling between the forward and 
backward propagating waves is proportional to d a \np n times an oscillating phase. Thus, 
naively, one would think that backscattering effects are of order d a lnp n . However when 
the effective potential V n (a) is slowly varying (i.e. in the WKB limit d a p n /p 2 n <C 1) 
this is incorrect because successively backscattered waves destructively interfere. This 
destruction is so effective that the backscattering amplitude is exponentially small. 

To establish this, suppose that there are two asymptotic regions, a < a_ and a > a + 
wherein V n = const and hence the WKB approximation is exact. Let us consider the 
solution which initially contains no backward propagating wave: T> n (a < a J) = and 
C n (a < a J) = 1. The backscattering amplitude is then V n (a > a + ). It can be calculated 
perturbatively by taking C n = 1 in eq. @ (for a more rigorous justification of this result 
see H @): 

B n = V n (oo) ~ J da I^ e +2*r M«0 ~ e 2i T d -' Pn(a') (10) 

The second approximation follows from the fact that in the WKB limit, the integral is 
dominated by the (complex) saddle point a* where p(a*) = 0. For instance when V n (a) 
has a minimum at a = a : V n (a) = v n + ^f-(a — ao) 2 + 0((a — ao) 3 ), the saddle is located 
at a* ~ a ± and the backscattering amplitude is D„(+oo) ~ e -^ v "/2G^/v^ U p 0n 



neglecting cubic and higher order terms in a — do- The lessons from this exercise are 
that: 

• The naive estimate of the backscattering amplitude is incorrect. A more detailed 
analysis shows that it is exponentially small. 

• If the kinetic energy of gravity (= V n (a)/G 2 ) has a minimum at a , then the 
backscattering arises from a region of width Aa = (^-5 2 p n ) -1 / 2 around the mini- 
mum. 

However, when the effective potential V n (a) does not have adequate asymptotic re- 
gions, it is impossible to isolate unambiguously a backscattering amplitude. Indeed, 
since the WKB approximation is no longer asymptotically exact, one needs an addi- 
tional principle to distinguish backward from forward classical motion. This situation 
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generally arises in cosmology. For instance, for a de Sitter universe, V(a) = —a 2 + Aa 4 
and backscattering around any a ^> A -1 / 2 far from the turning point results from inter- 
ference over a distance Aa ~ a. So it is impossible to isolate the backscattering around 
a from the effects due to the turning point at A -1 / 2 . A similar situation obtains for a 
matter dominated universe V n (a) = —a 2 + 2aGE n : one cannot isolate the backscatter- 
ing around a given a from the contribution of the origin a = and the turning point 
a = 2GE n . In these cases, effects at the turning points (for instance boundary conditions 
at a = or a = oo) might dominate any backscattering which occurs in intermediate 
regions. 



3 The Wheeler-DeWitt equation in minisuperspace 

We now generalize the previous analysis to the case when the matter Hamiltonian is not 
constant. Our aim is to find the equivalent of eqs. (|]) so as to be able to reveal the 
interplay between matter transitions and backscattering effects. We thus consider 

G 2 d 2 a + V(a) + 2aGH M {ir, <j), a)] |*(a)> = (11) 

where Hm is the matter Hamiltonian operator. It depends on the matter coordinates 0, 
their conjugate momenta it and a. Notice that the wave function |^/(a)) is expressed as 
a ket only for the matter degrees of freedom. 

Our analysis of the solutions of eq. ( |TT| ) is based on a double hypothesis. First 
we require that the universe be macroscopic. This condition is satisfied if the radius 
of the universe is larger then the Planck length and if the total matter energy (Hm) 
is larger than the Planck mass. This condition guarantees that the propagation of 
a is WKB, i.e. < p. Secondly, we require that no individual matter quantum 

dominates the kinetic energy of gravity. That is we do not consider the case of an 
inflaton field whose expectation value (0) is macroscopic. Instead we are considering 
usual field configurations such that ((f)) ~ but with (Hm) ^ and macroscopic. 

This second condition breaks the otherwise existing symmetry between the radius of 
the universe a and the matter degrees of freedom. It ensures that a is the heaviest and is 
therefore singled out to parameterize the transitions among neighbouring matter states. 

We therefore postulate a neet separation of length scales: the Compton wavelength 
1/p of the universe is much smaller than both distance over which the momentum of 
gravity changes (dhip/ da)~ l (i.e. approximately the Hubble radius) and than l/(E n — 
E m ), the time scale associated with a typical matter transition. We emphasize however 
that out treatment is exact, i.e. we do not neglect corrections, we simply postulate that 
they are small. 

Since the matter transitions are the lightest, they should be treated quantum me- 
chanically. This is implemented by making an adiabatic expansion for the matter states, 
following |14| \TE\ . On the other hand the radius of the universe is the heaviest degree of 



freedom, and satisfies the WKB condition, hence we shall make a WKB expansion for 
the gravitational waves. This double expansion is the basic tool we use to analyze eq. 



7 



The adiabatic expansion for matter is realized by making the instantaneous (i.e. at 
fixed a) diagonalization of the matter Hamiltonian 



H M {a)\^n{a)) = E n (a)\if) n (a)) 
(ij; m (a)\^ n (a)) = 5 m , n . 



(12) 



We combine this with the WKB expansion in the following way. Let p n be the classical 
momentum if matter has energy E n (a), 



G 2 p 2 n (a) + V(a) + 2aGE n (a) = . 



Then we generalize eqs. (|5|, |^) and decompose as 

M*)) = E 

n 

E ~iPn{a 



e ~ i J a p n (a')da' g +i J a p n (a')da! 

C„ (a) , h VJa)- 



2p n (a) 



daMa)) 



g— i J a p n {a')da' 

C n (a) ; T> n (a 



2p n (a) 

g+i J a p n (a')da' 



2p n (a) 



2p n (a) 



(13) 



(14) 



These equations univocally define the C n and V Tl coefficients and ensure that these 
coeficients are constant in the limit wherein both the WKB approximation and the 
adiabatic approximation are exact. Moreover, using these expressions, the conserved 
current yields the simple expression: 



(*\id a \V) =EKU< 



\V n {a) 



constant 



(16) 



Taking the derivative of eq. ([14]) and comparing with eq. ([150 yields a relation 
between C n , T> n and their derivatives which generalizes eq. (|8|). Then inserting eq. (ISP 
into the WDW equation (^1|) and using this relation gives 

d a C n = Y,(da^n)^^^r^-^ da 'C m 
9aPn f a p n daj) 



+ 



2p ri 



+ E (da^mm^^^riP^da'^ 
m ^ V PnPm 



(17) 



and the same equation with C n T> n , i <-» —i. These equations are equivalentto the orig- 
inal WDW equation. The essential advantage of this rewriting is that it neetly separates 
the backscattering effects encoded in the last two terms from the matter transitions in 
the forward sector which are described by the first term. As an important consequence, 
it enables the WKB approximation to be implemented without factorizing, as usualy 
done in former analysis, a gravitational wave common to all matter states, i.e. without 
being obliged to make the hypothesis that matter is in a tight wavepacket in energy. 
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4 Unitary Evolution 



When the WKB condition d a p/p 2 <C 1 is fulfilled, there is no backscattering and one can 
correctly neglect the coupling between the forward and backward propagating waves. 
One thus obtains two uncoupled evolutions for forward and backward propagating uni- 
verses. The equation governing the forward sector is 

d a C n = £ (dM^^p^jP^'^Crn . (18) 

This equation has two nice properties [lj]. First it describes the unitary evolution^ of 
matter as a function of a and secondly it includes gravitational backreaction effects at 
the quantum level, i.e. not only in the mean. 

Unitarity directly follows from the fact that the right hand side of eq. flT8] ) is antisym- 
metric. This can also be deduced from eq. (|T6| ) and the fact that the coefficients C n (a) 
and T> n (a) evolve independently. To see more explicitly how this Schrodingerian char- 
acter is obtained, we define two effective matter wavefunctions associated respectively 
with the expanding and contracting sectors 

p—if a pnda' r 

= J2 C n(^~ trPnda '\Ma)) 

n 

= E^(«) e+JrP " da >n(a)) . (19) 

n 

They obey Schrodinger type equations 

±id a \cf>% f ) = H eff {a)\<j)% f ) (20) 

where the Hermitian operator H e ff plays the role of Hamiltonian. It is defined by its 
instantaneous diagonalisation in the basis \if) n (a)): 

( x/Pn — \ZpTrt ) ^ 

H e ff{a) nm = 5 nm p n (a) + i(d a ip m \ip n ) v " v ^— . (21) 

^ v PnPm 

The peculiar form of the second term arises from the interplay between the r.h.s. of eq. 
([18]) and the derivatives acting on the \i/j n (a)) in the l.h.s of eq. (|20"])p| 

Thus, in the absence of backscattering, <fj£ff have the following properties: 

• ^tff ls decoupled from (p^ff- 

• (j) e ff are local in the sense that they depend only on \^{a)) and d a \^f(a)) at a. 



1 Some authors jHJ] have advocated that the WDW equation should be replaced by a first order 
equation in d a so as to eliminate backscattering and guarantee unitarity. In our setting, this would 
correspond to legislate that eq. (|l^) is the "correct" equation governing quantum cosmology. 

2 Equation (21) corrects an error in |t4) 
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(fifjj obey a linear equation, since they depend linearly on the original ket |\&(a)). 
Thus they still obey the superposition principle. 

Their norm is constant since H e ff{a) is hermitian: 

(0+//(aM//(a)> = E\(Ma)\<Ptff)\ 2 = £ \Ua)\ 2 = const . (22) 

n n 

If one works with wave packets tightly centered around a mean energy E, cpfjj obey 

the time dependent Schrodinger equation idt(f)fff = H m(o±)0^// where H M (a±{t)) 
is the matter hamiltonian of eqs. (|TT| , |i2|) with a±(t) describing respectively the 
classical expansion or contraction of the universe driven by E. To see this, as in 
eqs. (0, d), it suffices to develop the state dependent functions p n (a) around the 
mean value p(a). To first order in E n — E, the second term of eq. ( ^IP vanishes 
and the first one gives 

H eff (a) nm ~ 5 nm {{p{a) + d E p{a){E n - E)} . (23) 

Up to a term proportional to the identity (=6 nm {p(a) — dEp(a)E}) which plays 
no physical role, this diagonal matrix is identical to that defined by the matter 
Hamiltonian, see eq. (|T2"D, since the overall factor dEp(a) = ±da±/dt is precisely 



what is needed to convert d a in <9j in eq. (|20|) . 



• When one drops the tight wave packet approximation, H e ff, defined in eq. fl2~T|), 
includes gravitational backreaction effects to the matter propagation which depend 
on n. These are encoded in the non trivial dependence of p n (a) on E n and in the 
second term on the right hand side. Then, matter evolution must be parametrized 
by a since it is meaningless to call upon a mean time parameter. 

In conclusion, we claim that (pf^ defined in eq. ([19]) are the unique functions of \I/ 
which have these properties. The only ambiguity lies in the definition of the WKB 
approximation since one could modify p n in eq. ([19|) by local terms of order p n (l + 
0(dlnp n /da)) without modifying these properties. This ambiguity reflects the choice 
of normal ordering of the operator d 2 a in eq. ([TT|). However, these modifications are 
negligible in the WKB regime. 

Before discussing the role of backscattering effects, we wish to recall that eq. (|IBD . 
the rewriting of is a convenient expression when matter evolves quasi adiabatically, that is 
when the distance scale over which matter makes transitions, given by id ln(E n — E m ) / da) 
is large compared to the wavelength (E n — £ m ) _1 . In this case one can easily calculate 
the non-adiabatic transition amplitude A m ^ n from \ip m ) to \ip n ). Moreover it proceeds 
along lines very similar to the computation of the backscattering amplitude T> n in eq. 
(prp|) . Indeed, one first neglects transitions and sets C n = 5 mn for all values of a. Then 
one inserts this zero-th order solution in eq. (|l^) and integrates over a to obtain 

A m ^ n = C n (a = +oo) ~ f da{d a ^ n )^^e^ a ^- Pm )da 

V VnPm 

~ e i P da' (p n {a')-p m (a')) ^ ^24) 



-1 
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In the second line, as in eq. (|10D, we have evaluated the integral by saddle point. The 
saddle is now at the complex value a solution of p n (a) = p m (a). 

Since a(E n — E m ) is much smaller than Gp 2 n (by virtue of the hierarchy of length 
scales discussed above), one can expand p n — p m in powers of E n — E m . Keeping only 
the first term in this series, one finds that the transition is given by the usual expression 
controlled by the imaginary time to reach the saddle point times the difference of energy 
of the matter states, see [111] for further discussion of matter transitions in quantum 
cosmology. 



5 Consequences of backscattering 

To determine the range of validity of the truncated equation (|T8|), one must compute 
the importance of the effects induced by backscattering. These effects arise when one 
abandons the WKB approximation. In this respect, it should be stressed that there are 
two types of corrections to this approximation: local ones which can be evaluated by 



expanding the wavefunction as a series in h ||16|1 , and global ones which mix forward 
and backward waves. Only the latter lead to drastic modifications of matter evolution. 
Thus, one should differentiate "in the WKB approximation" from "in the absence of 
backscattering" . Our computation of backscattering amplitudes shall be based on eq. 
(|l7l) and shall proceed as in eq. (|24D and eq. (pi]). 



Our first conclusion is that these amplitudes are exponentially small, whether or not 
matter dynamics are neglected. The new feature with respect to Section |2] is that in 
eq. ( |TTD there are now coupling terms between C n and all V m . This leads to the fact 
that the relative amplitude to backscatter into different states is approximately thermal, 
with the inverse temperature given by the twice the imaginary time needed to reach 
the saddle point. These backscattering transitions modify in a fundamental way the 
structure of the evolution of the coeficients C n since the evolution no longer closes onto 
itself. This implies two unusual features: first, an initially purely forward wave packet 
will inevitably leak out into the backward sector and secondly, the knowledge of the 
initial values of the D n is necessary to determine the evolution of the C n . 

Let us first analyse the backscattering without change of matter state. The second 
term of eq. ( |TTD has exactly the same form as the term studied in section ^|. Thus it 
can be analysed in the same way: let us suppose that the effective potential V n (a) = 
V(a) + 2aGE n (a) has a minimum at a n and that we can approximate G 2 p\ = V n (a) ~ 
V n [a n ) + ^f-(a — a n ) 2 - Then, if the universe is purely forward propagating in state n for 
a « a n , that is C m (a <<c a n ) = 5 nm and T> m (a <K a n ) = 0, the amplitude to backscatter 
into matter state n is T> n (a y> a n ). As in section this amplitude can be perturbatively 
calculated by saddle point technique. The saddle is the solution of p n (a*) = and it is 

located at a* ~ a n + i^J 2V ^"^ -. Upon neglecting higher order derivatives of V n (a), the 
amplitude to backscatter is 

B n _> n = V n (a » a n ) ~ exp ( _ 7 ^Yhl^ = J = exp f-^- Pn (a n )lm(a* n )] . (25) 

G \ 2V n V J 
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The quantity which appears exponentiated is the ratio of the distance over which V n (a) 
changes (= Im(a*)) to the wavelength of the universe (= l/p n ). It is the ratio of 
the longest length scale in the universe to the shortest. Therefore, p n (a n )Im(a*) ^> 
(Pn(a) — p TO (a))Im(a), c.f. eq. (|24l) , since the second expression is controlled by the 
difference p n — p m . Thus B n ^ n is exponentially smaller then A n ^ m . It is this latter 
inequality which determines the validity of the first order equation (|18|). 

The backscattering with change of matter state is mediated by the last term in eq. 



(17). It can be evaluated in the same way and is given by 



da (d a i) m \^ n ) Pn Pm e +t f d -' (Pn(a')+ Pm (a')) 



a* 

~ gif da ' (Pn(a')+p m (a')) ^26) 

where a* is the saddle point solution oi p n (a*) + p m (a*) = 0. To evaluate this saddle 
point integral, we define the mean energy E(a) = En + Em and expand the integrand to 
first order in E n (a) — E m (a): p n {a') +Pm(a!) = 2p(a') + (E n — E m )dEp- This yields 

|5_ m | 2 ~ e- 4Im ( r * da ' M x e- 2 /o mT * (27) 

where T* = J a da' 3eP is the time to reach the saddle point in the mean geometry. 
Thus the relative distribution of the backscattered states is approximately thermal with 
the extremely low temperature given by (2ImT*) _1 . 

In the above we have neglected the fact that backscattering can take place in several 
steps. For instance one can first backscatter from matter state C m to matter state T> n i 
and then change matter state from V n > to T> n . Such multi step transitions can compete 
with the direct transition. Determining which channel is dominant is a difficult task (see 
for instance p0|). Nevertheless a rough calculation of these multi step backscatterings 
shows that the result based on the direct channel gives a correct estimate: the amplitude 
to backscatter is an exponentially small quantity governed by the total momentum p n 
and the backscattered states are approximately thermally distributed. 

Note that if the (unusual) sign of the kinetic term of gravity were positive instead 
of negative, then the opposite would be true and backscattering to states with higher 
matter energy would be favored. Mathematically this change comes about because OeP 
would have the opposite sign. 



6 Interpretation: introduction 

The literature on quantum cosmology abounds with interpretations of the wave function 
of the universe This reflects the contradictory views that are held on what the 
wave function should describe and what it should predict. However the question of 
interpretation is, to a large extent, determined by mathematical and physical consistency. 

In order to develop a coherent interpretation we shall first consider simple physi- 
cal systems -electronic non adiabaticity in atomic or molecular collisions and particle 
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creation in external fields- whose mathematical formulation is extremely close to the 
WDW equation in minisuperspace models. For these systems the interpretations, i.e. 
the procedures that must be used to extract predictions from the theory, are known. 
They can thus serve as guides to suggest the interpretation in quantum cosmology but 
also as laboratories to formulate the restricted set of questions that we, living in the 
universe, can ask in quantum cosmology. We then return to the WDW equation and 
compare different interpretations of \I/ that have been proposed in the literature. 

The first conclusion of this investigation is that, upon neglecting backscattering, the 
current interpretation of Vilenkin appears to be singled out as the only consistent one. 
Indeed, it leads to probability amplitudes which satisfy the usual properties of quantum 
amplitudes: the superposition principle and decoherence of remote configurations are 
both guaranteed. This is not the case if one adopts the conditional probability inter- 
pretation. This interpretation instead does make sense for systems on which localized 
external devices can interact. Therefore it is inoperative in quantum cosmology since no 
such external system can be introduced. 

The second point concerns the physical interpretation of backscattering for the two 
examples. In each case, an additional principle, also based on the possibility of coupling 
the system to an external device, is required to reach the correct interpretation. In cos- 
mology, no such principle is available and therefore the interpretation of backscattering 
remains ambiguous. For instance whether or not one should "third quantize" and at- 
tribute some statistic to ^ would modify the physical consequences of the backscattering 
amplitudes. 



7 Electronic non adiabaticity in atomic collisions 

Let us begin by considering the problem of two colliding atoms. After factoring the 
center of mass coordinate and the total angular momentum, the residual Hamiltonian 
looks like 

p2 

Hatom = ^ + V ( R ) + H el{R, Qi, Pi) (28) 

where R is the relative distance between the two nuclei, P the momentum conjugate 
to R, and qt,Pi the electronic coordinates. If we consider an energy eigenstate of this 
Hamiltonian 

y E (R)) = E\v E (R)) (29) 

it has exactly the same structure as the WDW equation in minisuperspace, see eq. (pp . 
Indeed in both cases there is a heavy degree of freedom (the scale factor a or the nuclear 
coordinate R) and light degrees of freedom (the matter or the electrons qi). Because 
of this similarity the same techniques can be applied to both equations. 

The simplest approximation consists in treating the coordinate R classically as a given 
function of time: R(t). Then the residual degrees of freedom, the electrons, propagate 



r) 2 
°R 

2M 



+ V(R)+H el ( qi , Pi ,R) 
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according to the time dependent equation id t \ip e i) — H e i(R(t),q,p)\ip e i) and the non- 
adiabatic dynamics are encoded in the transition amplitudes from one instantaneous 
eigenstate to another one (Landau-Zener effect). 

If these transitions are too energetic or the motion of R is governed by a too spread 
out wave packet, the background field approximation is no longer correct. Then one 
must solve eq. (|29|) and the techniques discussed in Section |3] should be brought to 
bear|nj. Upon neglecting backscattering effects but taking recoil effects into account, 
one obtains a second regime in which the dynamics is governed by the equivalent of eq. 
(|20|) : Schrodinger like equations for the electrons in terms of the amplitudes C n (R) and 
T> n (R) associated respectively with forward and backward WKB waves governing the 
inward and outward motion of R. In these equations, unitarity is respected for the C n 
and V n amplitudes separately and the role of time is played by the heavy coordinate 
R. When backscattering effects in R cannot be neglected, these effective Schrodinger 
equations are no longer valid and one obtains the third situation wherein one must solve 
the full equations for coupled C n ,V n system, the equivalent of eq. (|T?D. 

The sole difference between the atomic problem and the WDW equation in minisu- 
perspace is the sign of the kinetic energy of the heavy nuclei which is positive whereas 
that of a is negative. This difference only plays a role in the third situation wherein it af- 
fects the identification of in and out modes and the spectrum of backscattered states, see 
the remark at the end of Section |^. As far as the light degrees of freedom are concerned, 
the second regime in atomic physics and in cosmology are completely analogous. 

For this reason it is very instructive to review the interpretation of the solutions of 
eq. fl25p. Indeed, the aspects which are often glossed over in text books are precisely the 
ones needed for the interpretation of quantum cosmology. In particular, we shall discuss 
two compatible but nevertheless different schemes of interpretation. In the first, one 
restricts the analysis to asymptotic amplitudes to find the system in a given state. In 
the second scheme, it is the wave function itself which is interpreted at all times. Upon 
investigating the relationship between these schemes, the separation of the length scales 
governing electrons and nuclei will again play a crucial role. 

The first interpretation, proposed by Born|21|, concerns scattering amplitudes. In 



the present case as well, eq. fl29|) describes a scattering event and the main question 
one wants to answer (and to compare to experimental data) is how the amplitude of 
finding a certain final state depends on the initial state of the atom. Let us therefore 
suppose that the atoms are initially in their ground state and approaching each other, 
i.e. C n (R = +00) = 5 u q. According to Born, the probability for the atoms to end 
up in state n is \T> n (R = +00) | 2 . The mathematical basis for this identification is 

the conservation of the Wronskian / dqi^* E id^ e- Indeed since the nuclei repel each 
other the region R < is inaccessible and one must impose ^e{R = 0) = 0. Hence 
the Wronskian vanishes everywhere (since it vanishes at R = 0). Therefore if one 
expands ^ e in terms of forward and backward coefficients C n and T> n as in section |^, 
the vanishing of the Wronskian implies that J2n l^n(-R)| 2 = J2 n \^n(R)\ 2 - In particular if 
J2 n \Cn(R = +00) | 2 = 1, then J2n [Dn(R — +°°)| 2 = 1 which expresses the conservation 
of probability during the scattering process. 

In order to deal with dynamically induced backscattering (which occurs in quantum 
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cosmology) it is necessary to consider situations in which the wave function is not re- 
stricted to vanish at a certain place. For instance one can think of a single atom moving in 
one dimension in a static but inhomogenous electric field. The Schrodinger equation for 
the atom has the same form as eq. fl29|) , but now R is the position of the atom and the R 
dependence of V(R) +H e i(R, q,p) is due to the electric field. Consider an atom incoming 
from R = — oo in a given energy eigenstate n . This is implemented by requiring that the 
solution satisfy the following boundary conditions \C no (R = — oo)| 2 = <5 n ri0 . Since there 
is no atom incoming from R = +00, one should also impose T> n (R = +00) = for all 
n. By these boundary conditions one specifies what are the in-modes, i.e. modes which 
possess a well defined semi-classical behavior in the remote past. In this determination, 
there is an identification of the asymptotic waves carrying a unit Wronskian with the 
corresponding classical trajectories obtained by the stationary phase condition applied 
to wave packets or by identifying directly the wave length times % with the momentum 
of the particle. 

Upon solving the full (untruncated) Schrodinger equation one obtains both the 
asymptotic amplitudes for the atom to continue over the potential barrier, C n (R = +00), 
as well as those to be reflected by the potential T> n (R = —00). Unitarity is expressed by 
the equality of the ingoing current to the sum of the outgoing currents: 

/ dqi$%id R V E = 1 = \C no (R = -00)| 2 = J2 \ V n(R = -00)| 2 + J2 \ C n(R = +°°)\ 2 • 

n n 

(30) 

We emphasize that in Born's interpretation the identification of the quantities T> w C n (R = 
±00) as scattering amplitudes, and hence unitarity, follow from the internal structure 
of eq. (p9|) and the use of asymptotic waves normalized to unit Wronskian for all n. An 
external concept was only necessary to identify the asymptotic incoming and outgoing 
modes (ie. the left and right movers). 

The essential simplification of this interpretation is that one considers only asymp- 
totic states. This leads naturally to the question whether one can also give an inter- 
pretation to the electronic dynamics in the intermediate region. The coefficients C n (R) 
and T> n (R) are natural candidates for this interpretation since they generalize to finite 
R the coefficients that where the basis for the asymptotic interpretation. Moreover, this 
is also supported by the fact that, when backscattering can be neglected, they obey a 
Schrodinger like equation with a Hamiltonian that tends in the limit of infinitely heavy 
nuclei to the time dependent electronic Hamiltonian H el (R(t),q,p). 

To further investigate this question, we now turn to the second interpretation of the 
solution of eq. (^). In this interpretation -which for reason of simplicity is generally 
the first one presented in text books-, \^e(R,<i)\ 2 is interpreted as the probability for 
the atom to be at R in electronic configuration q. The basis for this interpretation is 
two fold. First if one considers time dependent solutions ^(t,R,q) of the Schrodinger 
equation, then there is a conserved current 

d t \*(t, R, q)\ 2 + d R J R + d q J q = (31) 

where 

j R = md R y , j q = **id q * . (32) 
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Hence it is consistent to interpret |\&(t, R, q)\ 2 and therefore |\& e(R, q) | 2 as probability 
densities, see e.g. [^]. The second basis, due to von Neumann |p3fl, is the following: 
if we couple the atom to an idealized measuring device (initially in the state X R = 0, 
X q = 0) through a coupling 5(t)\q)\R)Px q Px R (q\(R\ where Px R and Px q are the variables 
conjugate to X R and X q , the measuring device will record outcomes R, q (i.e. X R = R 
and X q = q) with probability \^f(t, R, q)\ 2 . 

An essential point to note about these two arguments is that they both necessitate 
the introduction of external concepts which are used in situ and no longer only asymptot- 
ically: in the first case, external time, and in the second, the local measuring device. For 
this reason the application of this second interpretation will be problematic in quantum 
cosmology since these external concepts cannot be invoked. 

It is now very instructive to show that this interpretation -when used with care- is 
consistent with the S matrix interpretation. In particular, we shall show when and why 
the second interpretation also implies that the coefficients C n (R) and T> n (R) should be 
interpreted as the amplitudes for the electrons to be in state n when the nucleus are 
propagating forward and backward at (better near) R. In this, the necessity of having 
well separated length scales is crucial. 

Consider the probability \^>e{R, q)\ 2 that the atom is at R and the electrons in state 
q. In terms of the coefficients C n (R) and T> n (R) it takes the form 

/ p i J R PndR —i J R PndR\ 

\*E(R,q)\ 2 = Y.\\Cn(R)^— + V n (R) e )(q\MR))\ 2 • (33) 



This quantity is wildly oscillating because of the interference terms between C n (R) and 
T> n (R). These interferences arises from the localized measuring device which has inter- 
acted strongly with both the electrons and the nuclei. Since we are interested only in the 
electrons, we would like the measuring device to interact weakly with the nuclei. Thus 
we consider a measuring device that interacts with the nucleus over a certain range 
of R, and with a certain momentum sensitivity. That is we consider the probability 
that the nucleus be in a wave packet <Pr ,p = e~( R ~^ / A e tRP ° / V2tt A where as an 
example we have taken a Gaussian packet. If the interval A satisfies both A ^> l/p n 
and A <C (<9ijlnC n ) _1 , the probability amplitude to find the nuclei in state (pn ,p is 
approximately 

- -( Pn (flo)-Po) 2 A 2 

dRp Ro>Po (R)*E(R,q) ~Y, C n(Ro)(q\MRo)) = (34) 

J n V A Pn 

where we have neglected the R dependence of C n (R) {q\ip n (R)) over the interval A, and 
dropped the exponentially small contribution of the T> n coefficients. Thus only the 
amplitudes with momentum in the interval Pq ± 1/A are selected. From this expression 
it follows that 

2 Pn (R ) m 
is the probability that the electrons are in state n if the nuclei are in the interval Rq± A. 
Recalling we are working at fixed total energy E and that p n (Ro) is proportional to the 
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speed of the nuclei, it follows that we have calculated is the probability for the electrons 
to be in state n multiplied by the time they spend in the interval R$ ± A. This confirms 
the interpretation of C n (R) as the probability^ amplitude for the electrons to be in state 
n. 

In this derivation the identification of C n (R) as electronic amplitudes is inevitably 



approximate. Indeed eq. fl34|) followed from the fact that the wavelength of the nucleus 
and that of the electrons are very different. The same separation of length scales was also 
the justification for neglecting the coupling between forward and backward propagating 
modes in eq. (|17D , thereby obtaining a Schrodinger equation for the forward propagating 
coefficients only. Thus, this separation is twice used in order to reach the interpretation 
of C n (R) as the amplitude for the light degrees of freedom to be in state n. 

We now turn to the consequences of backscattering. In this case, J2 n \Cn{R)\ 2 is no 
longer constant. How is this interpreted in the context of electronic non adiabaticity? 

To address this question, consider a generic solution of eq. ( p9|) in the context of the 
atom propagating in a background field so that the whole real axis — oo < R < +00 
is accessible. For such a solution, none of the coefficients C n (±oo), £> n (±oo) are zero. 
Using the equality 

£ |C n (+oo)| 2 + £ \V n (-oo)\ 2 = £ |C n (-oo)| 2 + £ |P n (+oo)| 2 (36) 

n n n n 

the probability to be in state n and forward propagating at R = +00 is 

m(+oc)) = - ,„ 1 ^ (37) 



when the WKB approximation is asymptotically exact. This probability is a highly non- 
local concept since it involves the coefficients C n and T> n at + and — infinity respectively. 
For this reason it is a relevant and useful concept only for an external observer who has 
access to both the forward and backward waves at R = ±00. However, there are many 
relevant quantities which are local in R. Moreover these quantities govern the physical 
outcomes when one is in the second regime wherein backscattering is neglected. In ad- 
dition light subsystems can presumably only ask questions concerning these quantities. 
An example of such a local quantity is the relative probability to be forward propagating 
in state n or state m at R — +00: 



P(C n (+oo)) \C n \ 



P(C m (+oo)) |C m | 2 



(38) 



In such local quantities the absolute normalization involving a mixture of C n and 
V n coefficients disappears. However they do not obey a closed equation since the C n 
are coupled to the V n . Thus, the effect of the backscattering is to modify the effective 
Schrodinger equation for the forward propagating coefficients C n (R) into a stochastic 
equation, where the stochasticity is given by the absence of knowledge a forward propa- 
gating subsystem has about the backward propagating part of the solution. At the end 

3 For the same reason in relativistic theory one must remove 1 / y2w from scattering matrix elements 
in order to get probability amplitudes. This is known as the reduction formula. 
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of the calculation he must average over the unknown coefficients D n . Happily, when the 
electronic length scales are well separated from the nuclear length scales, this stochas- 
ticity is exponentially small. 



8 Relativistic particle in an external field 

In order to further investigate the relationships between the conservation of the Wron- 
skian, the identification of the asymptotic solutions and the interpretation of the wave 
functions, another situation should be considered: the Klein-Gordon (KG) equation gov- 
erning the propagation of a relativistic particle in an external field. This example has 
in common with the WDW equation a Lorentzian signature, i.e. the quadratic form 
determining the kinetic energy term is not positive definite. 

This analogy has been emphasized by many authors, see e.g. |24j for a discussion at 
the classical level, and has been used to advocate the necessity of performing a "third 
quantization". Indeed, in the case of the relativistic particle, unitarity is implemented 
by carrying out a "second quantization" . The difference with the non relativistic case is 
that the Wronskian is no longer interpreted as a particle current density but as a charge 
density. This follows from a different identification of the in and out asymptotic modes 
which in turn dictates the new interpretation. 

Let us review these features by considering the propagation of a charged particle in 
an external electric field. We shall take for simplicity the electric field to be homogeneous 
and to point in the x direction. We impose that the electric field vanishes for t —>■ — oo 
and for t — > +00 which provides us with asymptotic regions in which the momentum is 
constant and therefore the modes are exactly WKB. The Klein Gordon equation in the 
gauge A t = 0, A x = f(t), A y = A z = takes the form 

d 2 - (d x + ieA x (t)) - d 2 - d 2 + m 2 ]^ = . (39) 

Writing \1/ = e t( - kxX+kyy+kzZ ^ 4> k (t) , one finds that <f) k obeys the equation 

'd 2 + V k (t)}<p k = (40) 

where V k (t) = (k x + eA x {t)f + k 2 + k\ + m 2 . 

This equation is identical to that studied in section [| and can be analyzed in the 
same way. To reach an interpretation one first needs to identify the in-mode describing 
a particle incoming from t = —00. This is achieved by decomposing <ft k as C k (t)xk{t) + 
U k (t)x k (t), where Xkif) is the WKB solution with unit positive Wronskian of eq. (pEPp. 
The in modes are then given by imposing the boundary condition \C k (t = — 00) | 2 = 
l,T> k (t = —00) = 0. (Notice that the condition on the coefficient T> k is applied on the 
same "side" of the potential in opposition to what is done in the non-relativistic case). 
As in the non-relativistic case, the potential V k (t) induces backscattering effects and 
both the coefficients C k and T* k are non vanishing at t = +00. The conservation of the 
Wronskian now reads 

1 = f d 3 xty*id t ty = \C k (t = -00) | 2 = \C k (t = +00) | 2 - \V k (t = +00) | 2 . (41) 
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How should one interpret these coefficients at t — +00, and in particular the fact that 
Ck{t = +00) > 1? This is the content of the Klein paradox that was discussed in the 
early days of relativistic quantum field theory |25| . 

Probably the simplest way to reach a physical understanding is to proceed as in the 
non-relativistic case by building wave packets (in k) and following the classical trajecto- 
ries they describe. One then finds that the wave packets made of C^Xk follow, both for 
early and late times, the trajectories of a positive charged particle. But the part of wave 
packets proportional to V^Xk follow trajectories of negatively charged particles. Thus 
the conservation of the Wronskian expresses charge conservation and the decomposition 
into Ck and T>k is a decomposition into positive and negative charged particles. Hence 
the extra charge carried by the Chit = +00) is compensated by the presence of anti- 
particles that reach t = +00: \T>k(t = +00) | 2 > 0. Another more mathematical way to 
reach this conclusion is to note that eq. (|39|) is invariant under gauge transformations, 
and that the Wronskian is the charge associated to the gauge field A^. 

Once this is realized one still faces the problem of extracting (probabilistic) predic- 
tions from the theory. This is done by second quantizing the coefficient Ck becomes 
an operator destroying a particle of momentum k. In this framework, one establishes 
that in addition to the "elastic" scattering of initial (anti) particles, one is describing pair 
creation. Even in the in-vacuum state, pairs of particles are spontaneously created with 
probability amplitudes governed by the backscattering amplitude Bk = T>k(t = +00) 
calculated in eq. (|i0|), which in this context are called Bogoljubov coefficients. If some 
particle is initially present, one finds that there is also induced emission. In this rea- 
soning the conclusion that the Wronskian is expressing charge conservation relied on an 
external input, namely that t in eq. (^0]) is an external time which allowed to follow 
asymptotic motion. Moreover, the prediction that many pairs can be produced can be 
checked experimentally. 

We now discuss how these new effects modify the dynamics of internal degrees of 
freedom. The most natural way to implement internal degrees of freedom is to go back to 
eq. ( |39|) and take the momenta along the x, y, z directions to be the additional degrees of 
freedom. The model becomes non trivial when the electric field is not homogeneous, so 

ti licit kg j 

are not constants of motion (they can be replaced by adiabatic constants 
of motion as in section [$]). The choice of how the field should be second quantized is 
dictated by the spin-statistics theorem ||26||, ie. \1/ is either a bosonic or a fermionic field. 
The important conclusion for us is that these different second quantization procedures 
are not equivalent, even if one restricts oneself to the dynamics of the internal degrees 
of freedom. 

To second quantize one first selects a complete set of modes with only positive fre- 
quency at t = —00, (j) l k n to which one associates a particle destruction operator a\ n . To 
the complex conjugate of these modes one associates an antiparticle destruction opera- 
tor b\ n . If we neglect backscattering then the number operators N part = J2k a T a T an d 
Nanti-part = J2k &jfc™ commute with the second quantized Hamiltonian. The Fock space 
then decomposes into non interacting sectors. The simplest sector is the vacuum |0). 
The first non trivial states are the one particle (or antiparticle) states J2k /(^) a fc 
These states are exact analog's of the forward propagating modes in the atomic example 
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of the previous section, and should be interpreted in the same way. The next states are 
the two particle states y f(k, k')a\ n a™, ^|0) where / is symmetric or antisymmetric 
according to the statistics of the field. Since / does not factorize into a product of a 
function of k and k', the interpretation of this state cannot reduce to that of one particle 
states. Moreover the statistics play a role in the determination of the final amplitudes 
since the presence of a particular quantum modifies the spectrum of transition ampli- 
tudes through either Bose enhancement or Pauli suppression. Whether or not this has 
an interpretation from the point of view of an internal degree of freedom is unclear. 

The situation is worse when one takes backscattering into account. In this case, 
the number operators no longer commute with the Hamiltonian because of the non- 
vanishing character of the Bogoljubov coefficients (calculated in eq. fl25|) and 
(|2lf)) relating in and out modes. This implies that in order to determine the amplitude 
to find a given state k in the one particle sector, it is necessary to know all the initial 
amplitudes (in Fock space) to find particles and antiparticles. Thus the amount of 
initial data needed in order to extract exact predictions from the theory in the one 
particle sector is (infinitely) larger than in the non relativistic case. For this reason 
third quantization is not a very attractive alternative since it diminishes the predictivity 
of quantum cosmology. Furthermore the additional ingredients which imposed second 
quantization, namely the presence of an absolute time and of a gauge field coupled to 
\1/ are not present in quantum cosmology. For these reasons it seems advisable, in the 
absence of additional strong arguments in favor of it, not to resort to third quantization 
in quantum cosmology. 

9 The conditional probability interpretation based 
on the norm of ^ 

In the light of the previous examples we now return to the interpretation of the WDW 
equation. 

For many authors, conservation of probability is a necessary requisite for quantum 
cosmology to possess a coherent interpretation. The most widely adopted scheme which 
implements from the start unitarity and conservation of probability is the conditional 
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probability interpretation |27j In this interpretation one identifies 

<^ (a ' 0) = T77777 !\, 7 (42) 



as the amplitude of probability to find matter in state <fi at radius a. By construction the 
norm square of <p e ff is constant, independently of the validity of the WKB approximation 
or any other condition. 

The main problem with this construction is that e //, contrarily to (fyftf defined in 
eq. fll9|) , does not obey a linear equation (since \I/ does). Therefore it does not satisfy 
the superposition principle. This can be ignored as long as one considers tight wave 
packets such that the factorization of a single gravitational wave, common to all matter 
states, offers a good approximation. However, once the spread in matter energy is large, 
the probabilities significantly vary even in the total absence of interactions. To see this 
suppose that the matter Hamiltonian has only two instantaneous eigen-states | V^o) and 
Then, in the WKB regime, the wave function of the universe is 

-i J a dapo{a) -i J a dapi(a) 

1*00) = ft 1 , r- M+Ci i IVi) • (43) 

Let us further suppose that it is impossible to make transitions from state to state 
1. Then the instantaneous states are truly stationary: d a \ipo) = d a \ipi) = and the 
coefficients Co and C\ are constant. Using the conditional probability interpretation, one 
then finds that the probability for finding matter in state |t/>o) is 

P(V> |a) = ^il . (44) 

ICoP + ^ldl 2 

It reduces to the usual expression Po(a) = |C | 2 only when p (a)/pi(a) = 1, i.e. when 
E = Ei. Therefore Po(a) varies once ^ E 1 through the ratio of the two gravita- 
tional momenta p and pi entangled with the two matter states by the constraint. The 
more remote these energy are, the bigger the effect since the associated momenta will 
be more different. Moreover, upon taking into account backscattering effects, the prob- 
ability receives in addition to the continuous component a rapidly oscillating component 
with frequency 2p whose origin are the interference between expanding and contracting 
solutions. 



4 Although 0] does not use the name "conditional probability interpretation" , it implicitly adopts 
it. Indeed in |3[ the decomposition of the wave function into a gravitational and a matter wave is 
uniquely fixed by imposing that the matter wave satisfies J d0x*(a, 4>)xs{a, 4>) — const. Hence Xs(aj 4>) 
is identical to <^ e //(a, <fi) defined in eq. (p2[). 

Another recent paper that claims that matter evolves unitarily in quantum cosmology 28 contains 
mathematical errors. Indeed ^p{h a ) in eq. (7) must carry an index n since the equation it obeys depends 
on n (see [ flO| for a discussion of this n dependence) . But this implies that the second term on the right 
hand side of eq. (10) is non hermitian. Furthermore the third term proportional to f2 is non hermitian 
(it would only be hermitian if integrated over h a ). Hence the conclusion of J28| that ^ n |c„| 2 is constant 
is incorrect. 
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In our opinion, these two properties are sufficient to reject the conditional probability 
interpretation. Before explaining our reasons, it is interesting to recall why a similar 
construction can be meaningful for a non-relativistic atom. In that case, the probability 
to find the atom in electronic state under the condition that it be at R is indeed 
given by eq. (fE^). What does give a meaning to this expression is the possibility (and 
the necessity) of coupling the atom to a localized external device so heavy that its own 
recoils be negligible. Only then does the factor po/Pi i n e Q- take meaning as the 
relative time each state interacted with the detector, and only then can the forward and 
backward waves give rise to the interferences at frequency 2po- But if one couples the 
atom to a measuring device that does not significantly disturb the motion of the nuclei 
(see section |7|), one obtains the current interpretation based on eq. (^) instead that 
based on the norrr|], see eq. (|33|). 

To present our arguments against eq. (|44]) we first recall the following point. In 
quantum mechanics, from the solutions of the linear equation determining evolution 
(which can be either the Schrodinger or the WDW equation) one obtains that very 
remote states neither interact nor interfere. From this mathematical property and the 
conventional interpretation of probability amplitudes, one then gets that the probabilities 
of such remote states do not influence each other. But this fundamental decorelation of 
remote configurations does not obtain in eq. (|44|). 

From an epistemological point of view, the fact that linearity of the WDW equa- 
tion does not reflect itself in a simple property of the physical probabilities \4> e ff\ 2 is 
disturbing. Indeed if linearity does not manifest itself in a simple way at the level of 
observation, why insist that we start with a linear equation? We might as well introduce 
non linearity at the level of the WDW equation directly. 

In our opinion, the only way to save this interpretation probably requires to first 
restrict the applicability of eq. (|42| ) to neighboring states only. But then its status would 
be of an approximate character rather than a starting point from which all solutions of 
the WDW equations can be interpreted. 



10 Vilenkin interpretation within the WKB regime 

Of the different interpretations of the wave function of the universe which have been 
proposed, the closest to the physical examples presented above is that of Vilenkin[|J, 



as generalized in |T3[]. Indeed the essential lesson of these examples is that they give 
physical substance to the coefficients C n and V n . This is not surprising since these 
coefficients have the following mathematical properties: 

• The conserved current is expressed most simply in terms of these coefficients upon 
working with modes normalized to unit Wronskian, 



5 Some authors have proposed to introduce an additional variable rj such that the modified WDW 
equation reads id v ^ = H ty, thereby recovering a situation similar to that of an atom in a time dependent 
state. It is then very tempting to interpret \^f(a, <fi, rj)\ 2 as the probability to find a and 4> at r/. However, 
since in cosmology one must restrict oneself to questions such that only little momentum is transferd 
to a, one recovers the current interpretation of eq. ([38]) (see section Q). Thus this latter interpretation 
naturally arises irrespectively of the use of an additional variable as time. 
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• they obey a linear first order equation in a, 

• if backscattering is neglected they evolve independently and unitarily. 

In the absence of backscattering, these properties designate the C n and V n as the prob- 
ability amplitudes for matter to be in state n in the forward or backward sector. We 
invite the reader to compare this reasoning with the original one of Max Born, see ET 



Very important is also the fact that in this interpretation, decoherence effects arise 
naturally and obey usual properties, in contradistinction to what occurs in the condi- 
tional probability interpretation. For instance, when two states or two groups of states 
have only vanishing matrix elements of H e ff{a), their evolution is completely indepen- 
dent. It is thus fully legitimate to interpret the wave function composed of their sum as 
representing two different universes rather that a quantum superposition of states living 
in the same universe. This decoherence is the most basic one (when compared to that 
induced by the environment or measuring devices) since it is solely determined by the 
nature of the full Hamiltonian, i.e. the WDW constraint. It therefore provides a natural 
explanation for the fact that one should not consider quantum superpositions of states 
giving rise to different semi-classical geometries (histories) since these states completely 
decouple. 

The neglection of quantum backscattering effects can also be conceived as a mani- 
festation of decoherence. Indeed, the appropriate character of the decomposition into 
forward and backward solutions (rather than in sines or cosines) follows from the effec- 
tive decoupling of these two sets. Thus, there is a hierarchy of decoherence associated 
with the hierarchy of length scales. Indeed, the decoupling of the forward and backward 
sectors is the most efficient since backscattering amplitudes are exponentially smaller 
than matter transition amplitudes. This leads to a clear separation of the solutions into 
two sets characterized by the sign of p n {a)- 

Thus the only aspect which differs from that presented in the physical examples of 
Sections |7| and |8| is the identification of these sets with forward and backward propa- 
gation in time. In the previous examples, it was necessary to use an external time to 
identify which asymptotic modes correspond to forward motion at t = — oo. (Recall 
that this identification was different for the non-relativistic and relativistic examples.) 
In cosmology, an external time does not exist and thus the eventual determination of 
forward motion in time, i.e. the arrow of time, should be intrinsic. Presumably, the 
answer is thermodynamic: if the matter is out of equilibrium, then the arrow of time is 
determined by the direction in which entropy is increasing. This point of view was re- 
cently advocated in the quantum cosmological context in p9| . If this is indeed the case, 
quantum cosmology is incomplete since it can only be interpreted for macroscopical uni- 
verses containing enough degrees of freedom so that thermodynamics applies. For such 
universes, the hierarchy of matter and gravity scales that we twice exploit will be well 
established. Therefore it appears that only macroscopic universes can be meaningfully 
investigated and interpreted. 
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11 The interpretation in the presence of backscat- 
tering 

As emphasized by Vilenkin, since the WKB approximation is not exact, "probability" 
and "unitarity" are inherently approximate concepts. Indeed neither the norm of 4>tn 
nor of (j)~ff are separately conserved in the presence of backscattering. 

However this is not an insurmountable problem for the following reason. Any subsys- 
tem of the universe, for instance the degrees of freedom within the cosmological horizon, 
is necessarily an open system and its dynamics are therefore necessarily non unitary. 
The coupling between forward and backward propagating universes will appear to this 
subsystem as an additional coupling to the environment. Furthermore this coupling gives 
rise to exponentially small effects, with appearing in the exponential approximately the 
ratio of the Hubble radius to the wavelength of the universe. Hence it is completely 
negligible compared to the other couplings to the environment at least for macroscopic 
universes. 

Having said this, one may nevertheless want to know what the possible effects of 
backscattering could be -however small they are-, and how they should be interpreted. 

A first point to be noticed is that in most realistic cosmologies, backscattering cannot 
be dissociated from effects at the origin of the universe. The reason being that the 
distance over which the backscattering occurs is comparable to the distance to 

the origin of the universe. For this reason backscattering and its interpretation may be 
intimately linked to the ultra violet structure of quantum gravity. If this is the case, the 
interpretation cannot be discussed in the context of a truncated WDW equation. 

Nevertheless one may imagine universes in which asymptotic regions exist. Then 
backscattering events can be localized and their amplitudes properly evaluated. In this 
case, one can inquire into their interpretation even in the context of simplified WDW 
equations. As discussed in sections and § there are at least two consistent interpre- 
tations of backscattering. In both cases the choice of interpretation was dictated by 
reinserting the system into a wider context. In the cosmological context we do not 
know into what (if any) wider context the WDW equation fits, and thus the precise way 
backscattering should be interpreted is unclear. 

Let us nevertheless briefly discuss the consequences of backscattering. For a subsys- 
tem of the universe which is forward propagating, the effect of backscattering is included 
by averaging over the (unknown) state of the rest of universe. In a first quantized con- 
text, the rest of the universe would simply be the backward propagating waves. In the 
second quantized context it would correspond to all the other sectors of the theory with 
different "universe and anti-universe numbers" . Thus the averaging, and hence the effect 
of the averaging, is different in both cases. 

Moreover, the interpretation of backscattered waves themselves is problematic. For 
instance what is the physical principle which would determine the arrow of time in the 
part of the wave function whose origin is backscattering ? To be concrete, consider a 
solution of the WDW equation which for a <^ a contains only C n coefficients, but for 
a ^> a contains both C n and V n coefficients. (a is the center of the region where 
the backscattering occurs). From section || we know that the T> n (a) coefficients will 
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ultimately be more or less thermally distributed with a very low temperature. In which 
direction does the arrow of time point for the D n coefficients? Is a ~ a the origin or 
the end of the universe? 

In conclusion the separation between forward and backward propagating universes, 
and hence the appearance of an effective Schrodinger equation for matter appears to be a 
more robust concept than previously thought. Indeed the corrections to the Schrodinger 
equation are exponentially smaller than any matter transitions. On the other hand the 
precise way the backscattering manifests itself does not seem to be constrained by the 
WDW equation -at least in the simple models we have considered in this paper. Hope- 
fully this interpretational problem can be resolved once we have deeper understanding 
of the structure of quantum gravity. The inclusion of anisotropics and inhomogeneities, 
and more importantly of the ultraviolet sector of the theory, could lead to a unique 
consistent interpretation. 
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